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Compensation of linear control systems by the classical 
methods cf root locus and frequency response is adequat 
when there is a single variable. However when there are two 
Or more variables compensation of physical systems for 


desired system dynamics becomes a complex task [1]. 


In modern control engineering systems characteristic 
eguations with two or more unknown parameters result. In 
1959, Mitrovic's Method [2]j} was developed to portray the 
relationship between the roots of the characteric equation 
and two ccefficients of the system. This nethod allowed the 
two lowest ordered coefficients of the Trackers ct. C 
equation to be varied. Siljak expanced this method to tae 
Coefficient Plane to ailow the variation of any two of the 


coefficients of the characteristic equation. For examcle: 


n 
2 j k 
(oj eo = BS + \ as (1-1) 
1 5 k 
k=0 
k#1i,j 
wheres on are real constants 


Penden ar real Varjables ms 
1 a 


Application of this method resulted in constant zeta, omega, 


and Sigma curves on the B -B_ plane specifying the rocts of 
i J 


equation (1-1) for any choice of € andw . 
ih 


In 1964, Siljak [3] extended Mitrovic's work so that 
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the two variable paraneterS may appear linearly in tke 
coefficients of the characteristic equation. consider the 


linear case of the characteristic equation: 


n 
kK 
f(s) = y. ee = Q (1-2) 
k=0 
where: a = bQ «+c + d 1-3 
oS aie 
and: and are real variable parameters 


Application of this method resulted in constant zeta, omega 
and sigma curves for any choice of @ and 2) on the a-@B plane 


specifying the roots of equation (1-2). 


Siljak's Parameter Plane Method prevides the designer 
Wath a Sitiple procedure for Hactoritrg a characteristic 
SiiicteOn wana Glsplaying the resuitS in a parameter plane 
die#granm. It allows the designer to obtain information about 
Systen stabilitv and the affect of parameter variations and 
adjustments on Stability. In the general case the Parameter 
Plane Method can be appiied to any control problem in which 
1t is necessary to determine how the variation of one or 


more parameters effect the root location of a systen. 


In 1967, Bowie [4] disceavered that a limitaticn existed 
in parameter plane theory. While attempting to solve a 
Sixth order characteristic eaquation, Bowie found that he 
could not find a compiete set of roots in an area orf the 
parameter plane that should have contained only real roots. 
BY substituting the values of qd and 8 into the 
Characteristic equation he found the rocts. 1s was 
discovered that a Set of compiex roots existed corresponding 


£0 a €-G®) pair in the left-half S~plane. Further 
Be 
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investigation showed that this ¢€-G) pair formed a straight 
a 


line in the parameter plane. This was contrary to the then 


existing parameter plane theory which stated that a €-wW 
n 


pair will oniy occupy a point in the parameter plane. 


Bowie called this line of constant ¢€ -constanty the 
n 


SINGULAR IINE as it was formed by a singuiar value of 


€¢ and W =<*the singular line added a new dimension to the 
n 


parameter plane because it allows the deSigner to hold one 
root location constant while moving others to achieve 


desired system characteristics by varying @ and @ . 


12 





A. THE PARAMETER PLANE 


The Parameter Plane Method was developed by Siljak [3] 
aS an extension to #itrovic's Coefficient Plane. The 
parameter plane 1S a mapping of points in the s-plane into 
the @ - @ plane. The following is the mathematical 


development of that transformation. 


Consider the characteristic equation: 


n 
f(s) : (2-1) 
SS} = as ZS 
k 
K=0 
where: a (k=Oe loosen) are real vartapies 
a2 7 
and: S = -ct+tjw = -lw tjo Uaic (2—2)) 
n n 


Examination or figure 2.1 shows 


(2-3) 


and Vi- 2 = Sing 


By substituting equation (2-3), equation (2-2) may ope 


expressed as: 
Ss =W (Cos6 + jSing) (2-4) 
n 


k k 
and s =W (Cosk@ + jSink@) (2-5) 
ti 


ye 





+Tm 


n Th 


HeGuicuz. | 
A POINT IN THE S-PLANE 





Since 


9 = cos 0) (2~6) 

it follows that 
Cosk@ = Cosk{[Cos(-¢) ] i297) 
and Sank* = Siank[ Cos (-7) ] © (2-8) 


Equation (2-7) can be expressed as the Chebyshev function of 


eae first kind, eS) and Similariy equation (2-8) can be 


expressed in terms of the Chebshev function of the second 


Kind, Sin Uf[-3$ ]. Substituting equations (2-7) ana (2-8) 
, k 


equation {2~5}) becanes 


k K , 
Seeeecomts (-4) + goin 0 (-¢) J (2-3) 
n k k 
k 
where: Demy = (=!) (E (C) 
K K 

(2-10) 

d U €} = 1 ae! G0 ¢ 

2D , ; = {(~1) , 3 


Substituting equation (2-10} into equation (2-9), then 
Supstituting equation (2-9) into eguation (2-1), equation 


(2-1) becomes 


Jn n : 
£ (s) » Ww (-1) "2 isi ay (-1) "yo gy = 0 (2-11) 
S)= a - ; + JjSin a - = - 
k n , 0) : ye k a ee 
, K=() ; k=0 
For any non~zero value of both the real and imaginary 


parts of eguation (2-11) must go to zero independently 


therefore, 
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kK 0 


3 ; Sta tar 
aay ) , 63) = 


The Chebyshev functions of the first and second kind 


C22) 


are given by the recursion formulae 


i = 2T = oF 
< (3) * (3) 1 (3) 


k ke 
(2-13) 
a ey 8) 
where: T, (y) =e uo) ae 
Be 3 ole 
Ga) aa 


To eliminate BOS) from equation (2-12) it is neted that 


T = JU - U 2-14 
, oS) , 3) p71 62 ( ) 
then equation (2-12) may be rewritten 
n n 
d( 1)“ Wt U (YT) ] >. ( Tae wu ) 0 (2-15a) 
= a is = ze 2 
k n k y on J 
k=0 k=0 . 
n 
k+1 k 
yen a Wu {3} = 0 (2-15b) 
k n k : 


Examination of equation (2-15a) shows that its first part is 
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a linearly dependent solution to equation (2-15b) for ail 


non-zero values of zeta therefore equation (2-15) becomes 


n 
1 : a 0 
— a — 
Gale cve Us) 
k=0 
(2-16) 
n 
) ( 1)" "5 (¥) 0 
—_ a — 
k n k 5 
k=0 
For the Linear case; 
a =Qb + Bc + a (2-17) 
kK k k k 
where: b, c , and d are real constants 
k k K 
and: Q-and G ace the variable parameters 


Bm stitution Of Equation (2-17) into equation (2-15) yields: 


BG + € + D = Q 
1 Ae 1 
(2~18) 
and BQ+c G 0) 6 
Z 2 2 
where: 
Ty n 
B 7 A Xe B Tae a (Cv) 
1 =) aa 2 “L,I ee 
k=0. k=0 
n n 
Cc yy me a E Ne ‘y (T) (2-19) 
= =— < = = es = 
1 a a 2 x eee: 
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n 


n 
‘3 k k k kK 
° cae k an oy 8 Dm d,' k n ey 


Application of Craner's Rule for the soluticn of 
Sigaultanous linear eguations to equation (2-18) yields the 


desired parameter plane solution equation: 


oP. ~ cA ee ~ BD, 

A REN ss EE gs 2 6S air aw ty 2-20 
BC -~ BC GC 
i 2 2 1 


Siljak [3] defined the parameter plane as @ rectangular 
plot with as the abscissa and as the ordinate. Equation 


(2-20) gives q and @ as a function of Wand}. By fixing W=uy 
n st 


and varying W over the range from zero to infinity in 
n 


equation (2-20), a curve. is produced in the A-@ plane for 
Pema xed evalua Oo. Zeta. This Curve specifies the % and@ 
pairs which will cause equation (2-1) to have a pair of 
BeOtrs Corresponding to the fixed value of zeta. Sa tale cay 


Dyeetaxingd G) = @)  jand varying zeta from -1 to +1 in 
n ni 


equation (2-20), a curve is produced in the @-@ plane for 


the fixed vaiue of W .. This curve specifies the q and § 
ni 


pairs which cause equation (2-1) to have a pair of roots 


corresponding to the fixed valne of gy) . Thus equation 
n 


iazUyealtows the Mapping of complex points in the s-plane 


to the parameter plane. 


fo transform the real axis points in the s-plane into 
the parameter plane, replaces S in eguation (2-1) by -¢ to 


obtain: 


18 





nN 


k 
£(s) = ». ea Nn ay = 0 (2-21) 
k=0 


Now equation (2-17) is substituted into eguation (2-21) to 


obtain: 


n 
k 
». ( oe Teco) = Q (2-77) 
k=0 
Ppt artying equation (2-22) obtains: 
Bip spe (cj Dit) = 0 (2-23) 
where: 
n 
k k 
B(o) = ». (-1) be 
k=0 
ni 
k kK 
C(a) = (=1) 4 (2-24) 
k=0 
n 
k k 
is ». (-1) tae 
k=0 


For a given value of o the functions B(c), C(c), and D(ce) 
are constants and equation (2-20) resuits in a straight line 
in the parameter plane with a locus of point corresponding 


to the real roots s=-o. 


A digital computer program PARAM A was written Dy 
Nutting {5} to aid in solving the characteristic equation 
and plotting the constant zeta, omega, found at the end of 


this thesis and a brief explanation is found in Appendix B. 
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Consider a system whose characteristic equation is the 


fourth order polynomial: 


4 3 2 
f(s)=s-+8.5s +(50+28)s +(12.5Q+258+42.5) s+508+25=0 2~25) 


Figure 2.2 shows the parameter plane for the following 


values of zeta, omeca, and Signa: 


OMEGA ZETA SIGMA 
3.0 0.00 O25 
4.Q O23 0 1.0 
Det G25 ee 
6.0 0.40 3.0 
7.0 0.45 4.0 
8.0 Oe.0 ent 
9.0 1200 - 6.0 


Examinaticn of figure 2.2 shows the region of stabiiity 
pounded by J=0.0, Y=t.0, and @=0.0. All values of K and @ 

inside this region will result ina stable system. The 
roots of the system can be determined by selecting a point 
in the plane and reading the corresponding values of zeta, 
omega, and sigma. Zeta and omega will form a complex of 


roots so that: 


“" 


s= -Sw + jw )-5 
nh n 


anasthe Sigia poitt will form the real roots of the form: 


Although it may not be immediatly obvious, there are 
two Sigma peints through each point in the region. Carerul 


interpolation or the plotting of additionaiines is required 
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to obtain the exact roots of the systen. But? “sthas: is 
generally unnecessary as the designer will interested in the 


general location of a specific pair of roots. 


Four points have been selected from figure 2.2. the 
position of these points was measured with an engineer's 
scale and the resuiting vaiue of and were substituted in 
Peuation (2-25). The roots of the equation were then solved 
yet he wuSe Of a polynomial rootfinder program on tne IBN 
360/67. Tabie TI-1 shows the coordinates of the points, the 
Mmeoalting reets and the cerresponding values ci zeta, omega, 
and sigma. As can be sten the parameter plane provides a 
Simple means of determining a root location while variable 


parameters are adjusted. 


-+.461+55.814 


-1.6154+55.783 


-2.278#456.075 | 





POTNTS FeneM FIGURE 2.2 
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@—-scale: & units/inch 
B-scale: 1 unit/inch 


po ofhjges V5 7 


PARAMETER PLANE EXAMPLE 
FOURTH ORDER POLYNOMIAL 


Ze 





B. THE SINGULAK LINE 


As stated in Chapter I Bowie [4] discovered tne 
existance of the singular line in the parameter plane. In 
investagating the cause of singular lines Bowie found that 
all previous work in the parameter plane was correct for as 
far aS jit went. PcmeCwdeeeiot, NOwWwever, dive a full 
explanation of the case when both the numerator and 
denominator of equation (2-20) formed singuiar matrices. 
Siljak [3] did address the fact that this condition couid 
exist, put the idea of a singular line in the parameter 


plane had not, up to that time, hean conceived. 


Bowie found that in order for singular lines to exist 


in the parameter plane the following conditions must be met: 


Bl = Feo = 0 
12 ZA 
Ge eece, D= 0 (2-26) 
1 2 24 
Bie Be =O 
Vn | 2 
where: 
n n 
B “4 no ‘ UJ ) B { 1) U 
— _ oO = = 6) 
1 K ns ee : k he 
K=0 k=Q 
n, n 
C = - | ec 0 ‘é =)" - | ce U 2-27) 
1 oo S od oe 2 eo) Se S ie 
k=Q0 k=G0 
n h 
a= ied U De aa U (¢) 
‘jee ie Pet peas) Oe 2 ving “kns k ’s 
~ k=0 k=0 
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mid: W and ¢ are the constantw-constanty values. 
ns s S 


eemations (2-18) are solved for in terms of B the 


lt is two linearly independent equations: 


D B D B 
1 1 2 2 
Ga- = -=-d=-—--d (2-28) 
c C C c 
1 1 ze Z 
Equation (2-28) is the singular iine solution. ie 


Mea straight line in the parameter plane with a slope of 


me and it intercepts the -axis at 7 a All points on 
1 


Mine represent a Singular value of zeta and a singular 
e of omega and therefore map into a Single pair of 


Mex points in the s-plane where: 


s=-$W + jw \/ 1-3. ome?) 
5 As ns S 


Tt is obvious that singular lines do not exist for all 


m Only for the -G) pairs which will satisfy equation 


acteristic equations and when they do exist they will 
S ns 

26). To aid in determining the existance of singular 
28s and plotting them in the parameter plane the computer 
| 


Joan SINGULAR LINE has been written. 


Briefly, the program uses the reai constants of the 

ficients of the characteristic equation and the desired 

1e aS to solve for value for any real positive values 
Me in equation (2-26). They and W values are then 
ns . Ss ns 

dd to generate equations (2-27) and these results are used 


test foc the existance of Singular lines. If singular 
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region there will be two real roots and a pair of complex 
roots corresponding to the constant zeta-constant omega 


vaiues of the singular lines. 


A pair of points nas been Selected along two of the 
Singular lines. Points A and B are on the line for J =0.7, 


&) =13.28 and point C and D are on the line for Sso om 
n 


G) =12.309. Table Ii-2 shows the coordinates, the resulting 
n 


roots and tke values of zeta and omega. For Table II-2 the 
coordinates of the point were measured with an endineer's 
scale and the values of alpha and heta were substituted into 
equation (2-30) and the roots were computed uSing a 


polynomial reotfinder subroutine on the IBM 360/67. 


fable Ii-3 is for the same points, the only difference 
being iat the coordinates were computed using the 
sSlope-intercept method utilizing the printed output from the 
program SINGULAR LINE. 


Table II-3 shows that this metnod of obtaining systen 
roots provides an accurate srethod for obtaining singular 
lines and that the iines a for constant zeta-constant omega 
values, while Table II-2 shows that reasonably accurate 


results can be cbhtained uSing a minimum of equipment. 
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@-scale: 5 units/inch 
B-scale: 1 unit/inch 


Figure 2.3 


SINGULAR LINES CN THE PARAMETER FLANE 
FOURTH ORDER POLYNOMIAL 
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POLWES@ EPROM FIGURE 2.3 


TABE. 1 i~3 








OMEGA | ZcTA 


POD T ALPHA Bova 


COMPLEX ROOTS 





-9.298459.486 .700 
A 

-9.2984£359.486 F100 
B 

-4..9204311.27 ~490 
c 

-4.92043511.27 400 
D 


POINTS FROM FIGURE 2.3 
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To fully understand the relationship of singular lines 
to the parameter plane it 1S necessary to examine what the 


Singular line 1s and what it represents. 


The sSingulac line exists aS a unique situation that 
occurs while mapoing the s-plane into the parameter piane. 
The necessary condition that must exist is that the 
determinants of the matrices which form equation (2-20) must 
be singular. When this occurs, a discontinuity appears in 
the mapping from the s-plane to the parameter plane. The 


result of this discontiniuty is the singular line. 


The discontinuity occurs ror the singular value of zeta 
and the singular value of omega, and it maps into the 
parameter plane as a straight line from plus infinity to 


minus infinity. The equation of this line is 


D 3 D 3 
ti 1 2 Zz 

S e Cc c 
1 1 2 Zi 


The slope of the line is ree and the bata intercept is 
4 5 


-D yc. 
1 he 


Untii this writing it was thought that is a value of 
alpha and beta from any point aiong the singular line was 
placed in the characteristic equation, the system would have 
a pair of complex roots corresponding the the value of the 
constant zeta-constant omega of the Singular line. however 


it is shown in the following examples that this resuit does 





mot alvays OCCUL. 


Consider the following aircraft pitch stabilization 


system [6]. 


243) Gh) 6,0 HH G6) —2 


D. 


H6) 


Figure 3.1 
Pekar erp tr Ch SPABILIZATION SYSTEM 
where: @ is the pitch rate command 
g is the pitch tate 
eee om pune ang he 


7 is the elevator angle 


and: G (s) 1S a pole-zere compensator 
G (Ss) is a hydraulic servo 
2 
G (3) is the airframe 
3 


H(s) iS a rate gyro 
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The transrer functions of the system components are: 


a (e) = 
for 2 
s+6 


eee 
G4 {S) = 
s + 54.85 + 3201.01 
Kaaleicaem, (5 ted <2) 
G 3 (S) = ; 
S{s + 7258S + 124.614) 
H(S) = s 


The characteristic equation for the uncompensated system is: 


5 uy 3 Z 
pues oe. 3Ss)6 Uf 6 374H1.01S «6+ «6(3109265 + 68732.6K)s 
+ (398891 + i98779K)s = 0 (7) 


Pocse system “15 Unstabbeszor K>3.27. It 1s desired to 
comvensate for K = 5 and to avoid the naturel frequencies o£ 


the airframe, gq) = 3.24, and the hydraulic servo, qu) = 7.049 
n n 


The characteristic equation of the compensated system is: 


6 5 - X 
peeto(o2. 38 +6) s f (3701.0 1% + 62.38Q)s 
5 
fe stoo2.> + 61732.8K + 3741.01) s fms oooo lat 19875 hd 


Z 
+ 31092.58) s + (198751FA+ 3988916)s = 0 (323) 


This characteristic equation with K = 5 was used with 
the program PARAM A to obtain the parameter plane curves, 
figure 3.2, and with the program SINGULAR LINE to determine 
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the existance of Singular line. Two singular lines were 
meana fon Cach Zeta value of 0.1, 0.2, 0.3, 0.4, and 0.5, 
and one singular line was found for each zeta value of 0.6, 
tee, 0.8, and 0.9. 


zeta = 0.6 was chosen for compenSation as it has an 


omega of 12.74 radians and wouid result in a tairly fast 


response without cauSing resonance at the natural 
Erequencies of the airframe and the servo. The singular 
line for zeta = 0.6 is the broken line on the parameter 


plane of figure 3.2. 


Three peints were selected on the singular line. The 
Mommmesrewere; A{1#.56,50), 8(21.13,150)8 and €(27.69,250}. 
These values of alpha and beta . where put intc the 
characteristic equation and the roots were solved for on tne 


IBM 360 with a volynonial rootfinder routine. 


There were two real roots and two pairs of complex 
roots for each of the points. \Contrary to tne supposition 
that all points along the Singular line will produce a pair 
of compiex roots corresponding to the value of the 
zeta-omega of that iine, in this case such roots were not 
found. The singular line with zeta = 0.6 and  %mega = 12.75% 
Should have according to prevolus concepts result in a 


complex rccot pair located at 
s = ~7.645 + 370.194 (3-4) 
This waS hot found to be the case. The compiex roots 
that were found did not agree with the singular line roots, 


but they did agree with the parameter plane roots. fnese 


results are as £ollows: 
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Sipocalte:s 20 units/inch 
B-scale: 50 units/inch 


Figure 3. 


AIRCRAFT PITCH SAPD E ZATION SYS@EN PARAMETER PLANE 
WITH ZETA=9.6~OMEGA=12.74 SINGULAR LINE 
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POIN? CONPLEX ROOTS ZETA OMEGA 
A -8.0834743.718 .908 8.697 
~2.512+}69.00 036 69.04 

B ~7.9444 46.987 751 10.58 
-16.664364.41 oi 66.53 

C -7.871448.012 F105 1s 2s 
-20. 754359. 36 . 330 62.87 


identical results were found when values were taken 
from other Singular lines in the same relative area of tne 
parameter plane. Ali the roots were found to agree With tine 


parameter plane, but not with the singular lines. 


The initial response to these results was that the 
Singular lines that were found were incorrect. However, 
after carefully checking the program siat were found yere 
MicerLect,. However, after carefully checking the program 
Sistem iAR LINE and hand caiculating two of tna singular 


lines, tne program was deemed to be correct. 


In an attempt to find a solution te this dilemma the 
Singular line exanples given by Bowie were examined. In all 
cases che Singular lines were Located outside the region of 


the normal complex root parameter plane curves. 


In the belief that the constant zeta curves in figure 
3.2 formed a loop as omega increased, values of alpha and 
beta believed to oe outside the region were tried. The 
Straight lines connesting the zeta curves in the figures 
Shown in this thesis are the resuit of the DRAW routine and 
are not part of the parameter plane. They result when 
intermediate points are not in the area of the plot. Again 
three points on the same Singular line were selected. These 
points were D(76.93, W000) , E(142.57,20900) and 
F(208.22,3000). These points are not located on figure 3.2, 


but the location was computed using the slope-intercept 
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method. These results are as follows: 


POINT COMPLEX ROOTS ZETA © OMEGA 

D -7.725459. 542 .629 12.28 * 
~23.314351.08 2415 56.15 

E - 7.688449. 856 2615 5G ae 
~23.474549.55 ~428 54.59 

P -7.675449. 966 .610 12.58 * 
~23.5045349.04 .432 54.33 


It is evident from the above that cne pair of ccmplex 
roots, indicated by the astrisk, 1S approaching the value of 
those whicn shculd result due to the singular line and that 
the value of zeta and omega for that pair of roots is 


relativeiy insensitive to large changes in alpha and beta. 


In light ofr these resultS an examination of the 
parameter piane in the region of these points was conducted. 
For reasons of clarity oniy three values. of zeta were 
plotted. These were 0.0, 0.6, and 1.0. The resuiting 
parameter ;lane 1s seen in figure 3.3. As can be seen the 
constant zeta lines do not loop around as prevoiusly 
enought, Higumare NdiLscontimuous and tie Singular line 


approaches asymptotically. 


What has happened is that as omega was increased on the 
constant zeta curve and approched the singular point the 
matrix which forms the denominator of eguation (2-20) 
becomes very much smaller than the numerator causing the 
SOuetitemezeta, dine to go to plus infinity at the singular 
point asymptotic to the singular line. Once the value of 
omega iS past the singular point it reenters the parameter 
plane at minus infinity, asymptotic to the singular line. 
Since there is only one singular point for this value of 
zeta the constant zeta line will continue through the 


parameter plane approaching infinity once more as omega 
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approaches infinity. 


The reason that this occurs is evident when the mapping 
of the constant zeta line form the s-plane to the parameter 
plane is observed. The constant zeta line in the s-plane 1s 


continuous fern W =9 to GQ) =, therefore the mapping of this 
n ni 


line into the parameter plane must also be continuous. 
However there is a discontinuity at the singular point where 
the constant zeta curve in the paramter plane goes to 
infinity for a finite value of omega. The singular line 
provides a continuous mapping of the constant zeta curve as 
it traverses from plus infinity to minus infinity. This 
line exists at infinity, but not in the parameter plane as 


meeis known. 


In an erfort to examine the change in alpha and beta as 
omega was increased a printout of PARAM A was made. However 
due to the iteration of omega in PARAM A it was not possible 
emobtdain SUELTIGCient data around the singular point. Tn 
PARAM A for an omega of about 12 radians the iterative step 
is about 0.9 radians. Therefore the section of PARAM A 
which computes the constant zeta curves was modified to vary 
ciuega over 1 radian in steps of 0.001 radians starting 0.5 
radians kelow the Singular point and increasing to 0.5 


radians above the singular point. 


In order to see how small increases in omega affect the 


values of alpha and beta some of the points are shown: 


OMEGA ALPHA BETA 
125 De 88.56 1297.87 
12.342 111.31 1640.01 
12.442 149.21 2220.85 
12.542 226.06 ERCE ee 
12.642 451.58 6826.15 
PES ay 144447, 2200380. 
12.743 -66473.4 -~1012572. 
12.842 -454. 96 -~6980.32 
12.942 -230.44 ~3564 249 
13.042 -154.32 -2405.44 
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a—-scale: 


100 units/inch 
B=-scale: 500 units/inch 


Figure 3.3 
FARAMETER PLANE FOR AIKCRAFT ae STABILIZATION SYSTEM 
ZETA = Q.6 





As can be seen for large changes in alpha and beta in the 
area Of the parameter plane where the singuiar line is 
asymptotic to the constant zeta curves results in small 
changes or omega. Since oneja 1S verv insensitive to 
changes on alpha and beta in this region the cencept of the 
Singular line is vaiiad from the standpoint of the engineer. 
memes possible £0 vary a parameter along the singular line 
in this regicn and expect to maintain a constant value of 


zeta and cnrega. 


In lidght of the unexpected results obtained in example 
Iii it was decided to carefuily reexamine example II that 
was given in the previous chapter. The parameter plane was 
expanded around the area of the origin where the Singular 
lines pass through the complex root portion of the parameter 


plane. 


Figure 3.4 shows this area of the parameter plane. The 
Semracant Zeta Curves ere for vaiues of 0.0, 0.2, 0.5, 9.6, 
(eo, and 1.0. BickecwdGeomtvOpcingudiar lines for each zeta 
value of 0.2, 0.5, 0.6, and 0.8 which are shown as dashed 
lines. There are two discontinuities in each zeta curve. 
Although it iS not immediately obvious due to the Limited 
area of the plot, all the constant zeta curves originate in 
the second quadrant, and become infinite at the first 
Preontruity end then réappeat at infinity in the fourth 
euadman t. The curves, for values orf zeta up to and 
including 0.5 then follow a path just to the left of the 
Origin and enter the first quadrant where they again become 
@iscontinuous and intinite. The zeta curves for values of 
zeta greater than 0.5 come from infinity in the fourth 
quardant and go ‘+o the left of the crigin where they turn 


mito the third quadrant where tney again become 
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Mmisconciauouwe amd iinzinite, The curves for zeta values of 
0.5 and lessS reappear at infinity in the third quadrant then 
[eee Oe ameintcy as a COntinuous function in the first 
quadrant. Tne curves fcc zeta values greater than 0.5 
reappear at infinity in the first Guadrant, Sweep toward the 


Geeayn then back through the £Z1rst Guaarant to infinit 
J y 


The region of Stability in thé parameter plane is 
bounded by the zeta=9%.0 curve, therefore the Singular lines 
of interest lie in the first quadrant. As can be seen the 
zeta=0.5 and zeta=0.6 Singular lines nearly coincide with 
the corresponding zeta curves in the parameter plane and tine 
zeta=0.2 and zeta=Q.8 lines are being asymptotically 


approached by the corresponding parameter plane zeta curves. 


AS in example III the modified section of PARAM A was 
used to abté@in a printout of alpha and heta for a 1 radian 
change in omega aS the constant zeta curve passed tarough 
the Singular line vaite. In this example the constant zeta 
curves in the first guadrant approeched their singular value 


very rapidly. 


For purposes of illustration points have been placed on 
the constant zeta curves in figure 3.4 at tne value at which 
omega iS within 0.5 radians of the Singular point which is 


Geeanrinaty. These points are: 


POINT ZETA ALPHA Bela 
A Ca20 4.078 5 ol 
B Oreo 0 =r 2 Uno 
Cc 0.60 0.577 1.60 
c 0.80 4.754 6.62 
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As can be readily seen these points are within 5 
percent cfr the Singular vaiue of omega and if one continues 
Pmetge tia Zeta Gumves iteill be feund, that for practrical 


design purposes, the changeS in omega are inconseguential. 


It has been shown that a true Singular line in the 
Original concept aces not exist. In examples III and IV it 
is shown that, i£ a designer carefully examines the Singular 
lines and their relationship to the constant zeta curves of 
the parameter piane the Singular concept can be used for 
practical design applications. 

EXAMPLE 


This example is the original system in which Bowie 
found that singular lines existed in the parameter flane. 
The system which was used has a characteristic equation 


which is the following sixth order polynomial: 


a 


-_ 


6 5 ; 4 3 Zz 
s + 80s + (2044+1600)s + 840ds + (1600844008) s 
+ 16006s + 16006 = 0 (325) 


fhe parameter plane with constant zeta curves with the 
memes ot 0-0, O-.2, O24, 6.6, 0.2, and 0.9 is seen in 
Piolirteswiaemand Se6. Figure 3.6 as the area anound the 
Setory an -fagure 3.5. There were t uo Singular lines tor 
PaermOumene geta values 0.2, 0.4, 0.6, 0.8, and 0.9 and 


these are shewn on figures 3.5 and 3.6 as dashed lines. 


As can be seen the parameter pliane has no visibie 
discontinuities as omega increases, however due to the fact 
that two singular lines exist for tach vaiue of zeta there 
must be two discontinuities in each zeta curve. Examination 


of the singuiar iines ana the constant zeta curves of the 
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parameter plane indicates that the Singular lines are 
murgent tc the Corresponding zeta curve at one point. when 
the values of zeta and omega for these points were examined 
it was found that the tangent points occur at the singtlar 


points in the constant zeta curves. 


Unlike examples III and IV where the mapping from the 
S-Dlane approached infinity as the singular pcint was 
approached this example maps as though no discontinuity 
existed. When the discontinuity is reached the constant 
zeta curve immediately goes to plus infinity tangent to the 
curve at the point of discontinuity and returns from minus 
infinity tangent to the constant zeta curve and continues 


Mapping to infinity as omega approaches infinity. 


This phemonenon does not reveal itseif in a normal 
parameter plane study for two reasons. The first of these 
is the previously mentioned size of iteration in the PARAM A 
program and the second reason is that a normal practice in 
colputer programing 1S to limit the numerator of a term to a 
minimum size to prevent an overflow or underflow in the 
computer when a number beyond the computer's capablity is 
reached. This practice was followed in the writing of PARAM 
A so if the Singular point was reached on an iterative’ step 
the program would have recognized this as a forbidden number 
and continued to the next value of omega, thereby bypassing 
ries singular pornt and coOntinting with no indication that 
this condition existed. If Bowie had not opeen interested in 
the area outside the normal complex root area of this 
example he would have had no indication that singular lines 


existed. 
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In this example the singular lines are true Singular 
lines as originally conceived. If a value of alpha and beta 
is selected at any point along these lines a pair of complex 
roots corresponding to the constant zeta~constant omega of 
that line will result. Yo show this fact six points have 
been selected in the region of Singular line on figure 3.5. 
As can be seen by the following the expected complex roots 


were obtained. 


SINGULAR POINT | COMPLEX ROOT 
POINT ZETA OUEGR EXPECTED OBTAINED 
A 0.80 Vib Pee he, —1G2¢445 14.13 —-18.424514.15 
0.40 “7.44 VO. 964 903.515 ~18.95435343.44 
B Cece 20200 = tome 2+79.101 Ae t a ocelot 
0.60 Sale oS Ss Ot 6 22 Sloe? Ite oS 
C 0.80 Sao oG Selig 2. 07S me ee OO 
Oo Sp wn 3/5) cs 2) (Jane A ear Oe 9 tt 42522 
D 0220 SUR sie ~3.450+350.728 “moeroot yO. 29 
0.60 es Seo | De 22 165 SUR Basle 745 
E 0.60 Vip oI Eero 2 UO leo a 0 2 
0.40 u7.47 =e ot 44 3.5 =e ott 743.43 
F 0.60 Jee | =e et 25 OS = (| Rag Janay ey) 724) 
0.30 Figo ae eo, = TO yd. 13 Saul ieegmicls so 


In summary, the singular lines result when a constant 
zeta or constant omega curve is mapped from the s-plane to 
the parameter plane and a discontinuity occurs. Ate thas 
discontinuity the parameter plane mapping goes to infinity, 
either plus or minus, and returns from the other side of 
Pie LnLt y . The singular line simply provides a mapping 
through this region to maintain a continuous mapping from 
the s-plane to the parameter plane. ieee se ss ven nas 
examples TiIT and IV, where the mapping first approaches 
infintiy as the singular point is approached the singular 
line concept iS a valid engineering aid over a limited range 
of the singular line. However in cases such as example V 


the singular line may be used over itS entire range as an 
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fee PARAMETER SPACE 


The parameter space is an extension of the parameter 
plane when there are three or more variable paraneters [7], 
moo) | 9 j- Once the concept of the parameter plane is 
understood it is relatively easy to extend the concept to 


three or more variable parameters. 


The easiest method to generate a three dimensional 
parameter space is to set the third variable parmeter the 
third variable parameter to a known value and then generate 
the two dimensional parameter plane. The process is then 
repeated for additional values of the third parameter. The 
result is a series .of parameter planes Which will form a 


MabkanmeteL Space. 


EXAMPLE VI 


Consider the following systen: 








K 


Ss (st 1) (S+2) 








Figure 4.1 
PieOeOnlGheoYot fii ALTH DUAL CASCADE COMPENSATION 
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iiwemenatacteristic equation for this system iz: 


5 4 s Z 
aeet (32 +G)s jemi + (60 + 926 4 Kis 


Pom AK + 50 )s + KA = 0 (4-1) 


Peeure 4.2 Shows the parameter space for the constant 
zeta=0.5 and constant omeya=71.7 and the values of K from 10 
Memioe It Was constructed by fixing the value of kK and 
computing the parameter plane uSing PaRAMN A. The curves were 
then plotted onto a single plot to show the three 


dimensions. 


Both the constant zeta curves and the constant omega 
curves frorcm a surface and the interconnection of these 
surfaces forms a curved line through the parameter space. 
Any point on this line will generate a pait of complex roots 
corresponding to the s-plane point of zeta=0.5 and 
omega=1.7. For these values of zeta and omega the conplex 


roots are: 


S$ = -0.850 + j1.470 (4-2) 


By selecting three points on this line this is shown to 
be true. MiemPOdiiis are A{O. 16, 2.6, 10.0), BOs ia a0 
ieee and €C (0.52, 3.63, 15.0).- Theses points (Were taken 


BEOte tie parameter Diane using an engineer's scale. 


POINT COUPLEX ROOTS SETA CREGA 
A —0.847+31.478 2497 1.703 
B ~0.849+57.475 498 1.701 
C -0.851451.469 2501 1.698 
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Pepe owJNGULAR SURFACES 


The parameter Space 1S COntinuous throughout the range 
of the third variable parameter. If for any value of the 
third parameter a discontinuity appears in the parameter 
plane it will also appear in the parameter Space as a 
discontinuity. AS stated in chapter [II this discontiauity 
Will result ina Singular line. As the third parameter is 
varied the discontinuity will result in singular lines for 
the same value of zeta and omega. The slope and the axis 
intercept points will change slighty so the result is a 
smooth surface representing a warped hypersurface traversing 


through the farameter space. 


Pere seoels wa protograph Of a Singular surface todel. 
The Singular surface is from the system in example VI. the 
Singular values are zeta=0.5 and omega=12.04 and the third 
Variable K 1s varied from 5.0 to 20.0. AS can be seen the 
surface 1S a smooth warped hyperplane. The model was 
constructed by obtaining the singular lines for the same 
Singular value at each value of Kk, transferring then to 
Garapoard, cutting then and accentting them with a 
contrasting color. The photography was done hy the Photo 
Division, Education Media Department, Naval Postgraduate 
Ssenool. 


Figure 4.4 is a drawing of the same singular surface as 
figure 4.3. The difference between the drawing and _ the 
photograph is that in the photograph the K axis starts at 
K=5 in the foreground continues back through to K=20, while 
in the drawing K=5 is the most distant singular line and the 
Singular lines progress forward until the value of K=20 1s 


reached. 
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The relationship of the Singular surface to tne 
parameter space is exactly the same as that of the singular 
line to the parameter plane. AS with the Singular line, the 
Singular surface provides continuity at infinity. The 
designer must be aware of the type of discontinuity which is 
producing the surface. Ii the discontinuities occur like 
the ones in example ILI and IV then the principfles .of 
Singular lines may only be applied in limited areas cf the 
surface. If the singularities occur like those in example 
IV then any point on the Singular surface will give a pair 
of complex roots corresponding to the constant zeta-constant 


omega of the surface. 
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Veeco urnNoe ALTON TECHNIOUES USING THE SINGULAR LINES 


omy eee ce es ee ee Se ee ee re gee ae ee ee ee — oe eee ee ee ——— os oe sew ee ae ee ee Se —s eo ee ee 


The major advantage of using Singular lines for 
compensation is that aS one moves up and down the singular 
line, Within previously mentioned limits, one pair of 
ccmplex roots remains in the same locations in the s~fplane. 
Wile tits palr Of FOOtS 1S stationary other root locations 
can be varied as the location on the Singular line is moved 
through other Singular lines, real root lines, and other 


constant zeta parameter plane curves. 


There many excellent textS available which go into 
great detail on how root location will affect systen 
response, overshoot, damping and so forth [1], [10]. Once 
the basic principles on how root locations affect a system 
are understood it 1S possible to use singular line concepts 


Momcotarn proper root Locations. 
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EXAMPLE Vill 


ee ee ae Eanes ee ae he ane 


Consider the pitch stabilization systen given in 


example I1il. 


ac S98 


- H6) 


Ergimes. 4 
PekG weer £ CH STABILIZATION SYSTEM 


The characteristic equation for the system is: 


6 5 uy 
Set to 2. 38 +3) s + (3741.01+ + 62.388)s 
3 
fms 2.5 +£.°61732.8K + 3741.018)s sees leet 6 7S 1 OX 
2 
+ 31092.56)Ss + (193751KKX+ 39838916)s = 0 (5-1) 


It waS shown in example IIIT that Singular iines exist for 
mbes yScen, oet trey are the result of discontinuities at 
Tee die y Therefore in order to effectively use the 
Singular line concert it is necessary to examwine the 
Singular lines to determine the areas in which Singular 


results could be expected. 


In the case of this problem it was decided to use the 


Singular line with zeta=0.6 ani omega=12.74 as a dominant 


2). 





root pair to give fast system response and avoid the natural 
frequencies of the servo and the airframe. Figure 5.2 shows 
the singular lines tor this system for values of zeta of 
ewe, 0.6, 0.7, 0.8, and 0.9. 


By moving up the zeta=0.6, omega=12.74 singular line 
above the region where beta=3000 the value of zeta and omega 
Should remain relatively constant. Two points have heen 
plmeected f£O attempt compensation, “A(273.87, 4000) and 
Peaouo.69, 48500). These values have heen put into the 
characteristic equation and the system was Simulated using a 


step input of 0.2 radians per second. 


Figure 5.3-is the response for the point A values of 
alpha and beta and figure 5.4 is the response for the point 
B values of alpha and beta. AS can be seen the systen 
response is essentially the same for both points. The 
response from the point A values has a peak overshoot of 
1.08 and the response from the point B values has a peak 
overshoot of 1.1. The steady state output is 1.66 degrees 


mmeeoth cases. 


This demonstrates one of the main advantages of the 
Singular line concept. As two parameters are varied the 


system response remains insensitive to the chance. 


If thas parameter variation were to take place inside 
the normal ccmplex root region of the parameter piane the 
Shift of root locations wouid result in a change of output 
response. Figure 5.5 is the parameter plane of the normal 
complex root area of the airdraft pitch stabilization 
system. Between point A and point C the linear change in 


ToewelomeuonamG the Change in alpha 1s 13.13. 


There are two pairs of complex roots and two real roots 


Wm this Section of the parameter plane. One of the real 
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memes 2s Cancemied at the origin of the s-plane by the zero 
produced oy the rate gyro in the feedback loop. fhe other 
real root location is pust s=-75 in these cases and will 
Only have a slight affect on the steady state value of the 
output. The complex root locations and their corresponding 


values of zeta and omega are: 


POINT CONLLEX ROOTS ZETA QUEGA 

A moeeoos jo. mig ~908 Seo 

ee oy © 9 0) BO) Sie: 69.04 

B sieht jo. 98/ Res fs Jn Be 
-16.66+7364.41 Sesh SOnoo 

C eo ia oO 12 wen Alis.2 3 

zn? Ot) 9.36 73 30 O2507 

HeOritietrosdonthant Foot 15 Located at s = -2.512 + 


Heemv0. At Point C the dominant root is located ats = 
oomcorvl + 36.012. ft can be seen that at the point A values 
the system response wlli be almost undamped with a peak 
mrerenoOt Cf about 1.9 and a natural frequency of about 69 
radians, while at the point C values the system will be well 


damped with a natural freguency of about 171 radians. 
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ale Corer STONSmwane RECONMENDATIONS FOR FURTHER STUDY 


PS a a we eS _—— a ae ea oe A A eee es a a a ee ee ee ee ee eee map EP coe wep RP eee ee ee a == ae 


It has been shown that some of the original concerts of 
the Singular line were in error due to previously unknown 
properties of the singular line and its relationship to the 
parameter plane. The main reason these properties were not 
found earlier was due tmainly to standard computer 


programming practices and limitations. 


Although it has been found that some limitations exist 
when using Singular line concepts it has been shown that the 
basic concepts of singular lines can be used by the design 
engineer to easily obtain desired system response when nore 


than one variable parameter is present in the systen. 


An example has been presented which shows that the 
Singular line concept is an excelient method to cbtain 
desired system response. The method will allow the designer 
to select either dominant root locations or, if it is 
desired, nondominant root locations. the root locations 
when selected from the Singular line region of the parameter 
plane will remain fixed in their location as long as the 
alpha and beta parameters are restricted to the singular 


Hane. 


FOr further study it is highly recommended that the 
case where true singular lines exist such aS in example V be 
thoughly studied to obtain a full understanding of why the 
Prt aemrotrmit Goes NOE OCCUr at infinity in the parameter 
plane. Once this is acconplished devise a method to force 
this type of condition at a desired zeta-omega value. liga 


this can be accomplished it would allow the designer to 
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emeerea FOOt Location in the sS-plane and use Singular line 


concepts for easy design and compensation. 


OZ 





— aS ee ee 


PARAM A 


This is a modified version of the PARAM A_ progran 
Written by Nutting [5]. Originally written in FORTRAN for 
the CDC 1604 Computer, it has been modified to FORTRAN IV 
Pores the L£BH 36076/ at Naval Postgraduate School. It has 
also been modified to allow the user to select the size of 
the output flot and to plot singular lines in the parameter 


plane. 


However PARAM A will not determine the existance of 
Singular lines. It is necessary for the user to know the 


value of ¢€and w of a Singular line before this opticn can 
n 


be used. These values can best be determined by use of the 
program SIfNGULAR LINE. The user then can either put these 
values into PARAM A or use the values to draw the singular 
line on the PARAM A output plot. 
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SINGULAR LINE 


The purpose of the program SINGULAR LINE is @) 
determine the existance of singular lines in the parameter 
Peemes Lhe program iS written in FORTRAN IV for the IBM 
360/67 computer and the CALCOMP plotter at Naval Postgrauate 


Senool. 


The values read into the program are the coefficients 
of the characteristic equation, the values of zeta to be 
examined, and plotting information. The program contains 
three subroutines which are located in the source library of 
the IBM 360/67. These are POLRT which has the ability to 
solve up to a 36 order polynomial, DRAW which outputs the 
meoce tO the CALCOMP plotter, and BCNV which converts numeric 


Symbols to alpha-numeric for use in DRAW. 


The program takes the values of zeta and first 
generates the Chebyshev function of the second kind. sD 
then generates the coefficients for the omega terms in the 


pcolynomial expansion of ae - pC = 0. This generation is 


acconplished by using the coefficients of the alpha and beta 
terms in the characteristic equation, the value of zeta and 
the Chebyshev function. The resulting omega coefficients 


are now input into the subroutine POLRT. 


POLRT solves for the roots of the omega polynomial and 
Gtepucasetne results to the progran. These roots are now 
examined for positive reali values. If positive real roots 


are present there iS a possibility of Singular lines 
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existing for that value of zeta. If there are not positive 
real roots, there are no singular lines. This information 


is printed out and the next value of zeta is examined. 


For the positive real roots Bre ear oo oe Dae and D 


of equation (2-27) are formed and used to insure that a 
Singular line does exist. If a Singular line does exist 


then Bite oa and a are used to generate the points along 


the line by using equation (2-28). These pointS are now 


input into the subroutine DRAW and are plotted and labeled. 


The printed output includes the value of zeta and 
omega, the alpha and beta intercept points, the slope, and 
the label. If there are less than two points generated 
Within the area of the plot there will be no plot of that 
imme. Chis information 1S indicated on the printed output 
and the user can either change the plot scales or use the 


Greput Jine Information to construct a singular line plot. 


This procedure is repeated for each value of Zeta, 
feertraing ali the tines on a Single piot. Multiple runs for 
different characteristic equations are possible by uSing 
additional data decks. Each equation will have its own 


Outreut plot. 
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